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We study cosmological application of interacting holographic energy density in the frame- 
work of Brans-Dicke cosmology. We obtain the equation of state and the deceleration pa- 
Q\ ' ramctcr of the holographic dark energy in a non-flat universe. As system's IR cutoff we 
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choose the radius of the event horizon measured on the sphere of the horizon, defined as 
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L = ar(t). We find that the combination of Brans-Dicke field and holographic dark energy 
can accommodate Wjj = — 1 crossing for the equation of state of noninteracting holographic 
dark energy. When an interaction between dark energy and dark matter is taken into ac- 
count, the transition of wd to phantom regime can be more easily accounted for than when 
resort to the Einstein field equations is made. 

. I. INTRODUCTION 

■ Recent data from type la supernova, cosmic microwave background (CMB) radiation, and other 

oc ; 

. cosmological observations suggest that our universe is currently experiencing a phase of accelerated 

>/~) ' expansion and nearly three quarters of the universe consists of dark energy with negative pressure 
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Nevertheless, the nature of such a dark energy is si ill the source of much debate. Despite 
the theoretical difficulties in understanding dark energy, independent observational evidence for its 
existence is impressively robust. Explanations have been sought within a wide range of physical 
phenomena, including a cosmological constant, exotic fields, a new form of the gravitational equa- 
tion, new geometric structures of spacetime, etc, see [2[ for a recent review. One of the dramatic 
candidate for dark energy, that arose a lot of enthusiasm recently, is the so-called "Holographic 
Dark Energy" (HDE) proposal. This model is based on the holographic principle which states that 
the number of degrees of freedom of a physical system should scale with its bounding area rather 

n n 

than with its volume [3( and it should be constrained by an infrared cutoff On these basis, 
Li [|| suggested the following constraint on its energy density p£> < 3c 2 m 2 /L 2 , the equality sign 
holding only when the holographic bound is saturated. In this expression c 2 is a dimensionless 
constant, L denotes the IR cutoff radius and m 2 = (87rG) _1 stands for the reduced Plank mass. 



Based on cosmological state of holographic principle, proposed by Fischler and Susskind 



6], the 
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HDE models have been proposed and studied widely in the literature 
HDE model has also been tested and constrained by various astronomical observations 
well as by the anthropic principle [3] . It is fair to claim that simplicity and reasonability of HDE 
model provides more reliable frame to investigate the problem of dark energy rather than other 
models proposed in the literature. For example, the coincidence problem can be easily solved in 
some models of HDE based on the fundamental assumption that matter and HDE do not conserve 

n 

separately [161 ] . 

On the other side, scalar-tensor theories of gravity have been widely applied in cosmology 
Q. Sea— tlMOri e S are not new and have a long history. The pioneering study on scalar- 
tensor theories was done by Brans and Dicke several decades ago who sought to incorporate Mach's 
principle into gravity [18]. In recent years this theory got a new impetus as it arises naturally as the 
low energy limit of many theories of quantum gravity such as superstring theory or Kaluza-Klein 
theory. Because the holographic energy density belongs to a dynamical cosmological constant, we 
need a dynamical frame to accommodate it instead of general relativity. Therefore it is worthwhile 
to investigate the HDE model in the framework of the Brans-Dicke theory. The studies on the 
HDE model in the framework of Brans-Dicke cosmology have been carried out in 19|, [2CJ, |2l| . The 
purpose of the present paper is to construct a cosmological model of late acceleration based on the 
Brans-Dicke theory of gravity and on the assumption that the pressureless dark matter and HDE 
do not conserve separately but interact with each other. Given the unknown nature of both dark 
matter and dark energy there is nothing in principle against their mutual interaction and it seems 



very special that these two major components in the universe are entire 
this possibility is receiving growing attention in the literature 
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y independent. Indeed, 



261 ] and appears to 



be compatible with SNIa and CMB data [27(. On the other hand, although it is believed that 
our universe is spatially flat, a contribution to the Friedmann equation from spatial curvature is 
still possible if the number of e- foldings is not very large Q]. Besides, some experimental data 
has implied that our universe is not a perfectly flat universe and recent papers have favored the 
universe with spatial curvature 28]. 



In the light of all mentioned above, it becomes obvious that the investigation on the interacting 
HED in the framework of non-flat Brans-Dicke cosmology is well motivated. We will show that 
the equation of state of dark energy can accommodate wd = —1 crossing. As systems's IR cutoff 
we shall choose the radius of the event horizon measured on the sphere of the horizon, defined as 
L = ar(t). Our work differs from that of Ref. |.19j in that we take L = ar(t) as the IR cutoff not the 
Hubble radius L = iT -1 . It also differs from that of Ref. [2oj], in that we assume the pressureless 
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dark matter and HDE do not conserve separately but interact with each other, while the author 
of [3] assumes that the dark components do not interact with each other. 

This paper is outlined as follows: In section [Til we consider noninteracting HDE model in the 
framework of Brans-Dicke cosmology in a non-flat universe. In section HTT1 we extend our study to 
the case where there is an interaction term between dark energy and dark matter. We summarize 
our results in section IIVL 



II. HDE IN BRANSE-DICKE COSMOLOGY 



The action of Brans-Dicke theory is given by 

S = f d 4 x^ (-<pR + ^ v d^d v ip + L M J ■ (1) 

The above action can be transformed into the standard canonical form by re-defining the scalar 
field ip and introducing a new field cf), in such a way that 

*=k- (2) 



(3) 



Therefore, in the canonical form, the action of Brans-Dicke theory can be written [2 



s = J d 4 x^ (y-^4?R + ^d^dA + L M J , 



where R is the scalar curvature and <f> is the Brans-Dicke scalar field. The non-minimal coupling 
term cj) 2 R replaces with the Einstein-Hilbert term R/G in such a way that G~g = 27r0 2 /a;, where 
G e f[ is the effective gravitational constant as long as the dynamical scalar field 4> varies slowly. 
The signs of the non-minimal coupling term and the kinetic energy term are properly adopted to 

(-| ) metric signature. The HDE model will be accommodated in the non-flat Friedmann- 

Robertson- Walker (FRW) universe which is described by the line element 

ds 2 = dt 2 _ a 2 {t) ( + r 2 dn 2\ (4) 

\ 1 — kr z J 

where a(t) is the scale factor, and k is the curvature parameter with k = —1,0,1 corresponding 
to open, flat, and closed universes, respectively. A closed universe with a small positive curvature 
(Qfc ~ 0.01) is compatible with observations [28|]. Varying action ([3]) with respect to metric ((H) for 
a universe filled with dust and HDE yields the following field equations 

M fl2 4) -^ 2+ !^*= (5) 
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jV ( £ + a 2 + \) - - W + -) <P 2 = pd, (6) 

4:L0 \ a a 1 ) uj 2lo 2 \ lo J 

4, + 3H<j>-^-(*+H 2 + 4V = °< ( 7 ) 
2cj \a a 2 / 

where H = a/a is the Hubble parameter, /?£) and pr> are, respectively, the energy density and 
pressure of dark energy. We further assume the energy density of pressureless matter can be 
separated as pu = Pbm + Pdm, where pbm and pdm are the energy density of baryonic and dark 
matter, respectively. We also assume the holographic energy density has the following form 

PD = 4^' (8) 
where (f> 2 = uj/{2ttG c q). In the limit of Einstein gravity where G e fr — > G, the above expression 
reduces to the holographic energy density in standard cosmology 

3c 2 3c 2 m 2 

The radius L is defined as 

L = ar{t), (10) 
where the function r(t) can be obtained from the following relation 

r® dr = r°°dt = Rh 

Jo VI - fer 2 J a a ' 
It is important to note that in the non-flat universe the characteristic length which plays the role 
of the IR-cutoff is the radius L of the event horizon measured on the sphere of the horizon and not 
the radial size Rh of the horizon. Solving the above equation for the general case of the non-flat 
FRW universe, we have 

r(t) = -^=smy, (12) 

where y = VkRh/a. Now we define the critical energy density, p cr , and the energy density of the 
curvature, pk, as 

30 2 # 2 3fc<A 2 
P - = -to"' Pk = 4W • (13) 

We also introduce, as usual, the fractional energy densities such as 

_ p M ^PM n .* 

n * = ^ = ^?« (15) 
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For latter convenience we rewrite Eq. (|16p in the form 

HL = ir D (17) 

Taking derivative with respect to the cosmic time t from Eq. ()10|) and using Eqs. (|12j) and f) 171) 
we obtain 

L = HL + ar(t) = —^= — cos y. (18) 

Consider the FRW universe filled with dark energy and pressureless matter which evolves according 
to their conservation laws 

p D + 3Hp D (l+w D ) = 0, (19) 
P M + 3Hp M = 0, (20) 

where wd is the equation of state parameter of dark energy. We shall assume that Brans-Dicke 
field can be described as a power law of the scale factor, <p oc a a . A case of particular interest 



is that when a is small whereas u is high so that the product auj results of order unity 19]. 

This is interesting because local astronomical experiments set a very high lower bound on u 301 ]: 

in particular, the Cassini experiment implies that to > 10 4 [31I l^ . Taking the derivative with 
respect to time of relation <j) cx a a we get 

4> = aH<f>, (21) 

4> = a 2 H 2 (p + a(f)H. (22) 

Taking the derivative of Eq. (j5J) with respect to time and using Eqs. (fTBj) and (f2T|) we reach 

p D = 2H PD (a-l + cos y) . (23) 



Inserting this equation in conservation law (|19j) . we obtain the equation of state parameter 

1 2a 2^^ 
WD = --- Y -—cosy. (24) 

It is important to note that in the limiting case a = (u — » 00), the Brans-Dicke scalar field 
jecomes trivial and Eq. (|24p reduces to its respective expression in non-flat standard cosmology 



7] 



1 2^ 

«;£> = -- ^— cos y. (25) 
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We will see that the combination of the Brans-Dicke field and HDE brings rich physics. For a > 0, 
wd is bounded from below by 

1 2a 2 v / fT^ 



w D 



(26) 



3 3 3c 

If we take £Id = 0.73 for the present time and choosing c = 1 

Wn = — % — 0.9. Thus for a = 0.15 we have wd = —1. The cases with a > 0.15 and a < 0.15 



331 ] . the lower bound becomes 



should be considered separately. In the first case where a > 0.15 we have wd < — 1. This is an 
interesting result and shows that, theoretically, the combination of Brans-Dicke scalar field and 
HDE can accommodate = -1 crossing for the equation of state of dark energy. Therefore one 
can generate phantom-like equation of state from a noninteracting HDE model in the Brans-Dicke 
cosmology framework. This is in contrast to the general relativity where the equation of state of a 



noninteracting HDE cannot cross the phantom divide 



a. 



In the second case where < a < 0.15 we 



have —1 < wd < —0.9. Since a ~ 1/u and for uj > 10 4 the Brans-Dicke theory is consistent with 



solar system observations 3l|, thus practically a ~ 10 -4 is compatible with recent cosmological 
observations which implies wd — —0.903 for the present time in this model. In both cases discussed 
above wd < —1/3 and the universe undergoing a phase of accelerated expansion. It is worthwhile 
to note that since a ~ I /to and to > 10 4 , therefore for all practical purposes Brans-Dicke theory 
reduces to Einstein gravity as one can see from the above discussion. 
For completeness, we give the deceleration parameter 

a H 



(27) 



which combined with the Hubble parameter and the dimensionless density parameters form a set 
of useful parameters for the description of the astrophysical observations. Dividing Eq. © by H 2 , 
and using Eqs. ©, HZD, (HQ and ([22]), we find 

1 



* 2a + 2 
Substituting wd from Eq. ([23| . we get 
1 



[(2a + l) 2 + 2a(aw -1) + Q h + 3Q D w D } 



[2a + l) 2 + 2a(aw - 1) + tt k - (2a + l)Oo - -Q 3 ^ 2 cos y 



(28) 



(29) 



2a + 2 _ 

If we take O^j = 0.73 and ~ 0.01 for the present time and choosing c = 1, auj ~ 1, u = 10 4 and 
cos y ~ 1, we obtain q = —0.48 for the present value of the deceleration parameter which is in good 
agreement with recent observational results 



3411 . When a — > 0, Eq. ([29]) restores the deceleration 



parameter for HDE model in Einstein gravity 



25] 



[Id 
2 



3/2 
D 



cos y. 



(30) 
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III. INTERACTING HDE IN BRANSE-DICKE COSMOLOGY 

In this section we extend our study to the case where both dark components- the pressureless 
dark matter and the HDE- do not conserve separately but interact with each other. Although 
at this point the interaction may look purely phenomenological but different Lagrangians have 
been proposed in support of it [35]. Besides, in the absence of a symmetry that forbids the 
interaction there is nothing, in principle, against it. Further, the interacting dark energy has been 
investigated at one quantum loop with the result that the coupling leaves the dark energy potential 
stable if the former is of exponential type but it renders it unstable otherwise [36]. Therefore, 
microphysics seems to allow enough room for the coupling. With the interaction between the two 
dark constituents of the universe, we explore the evolution of the universe. The total energy density 
satisfies a conservation law 

p + 3H(p + p) = 0. (31) 

where p = pm + Pd and p = pd- However, since we consider the interaction between dark energy 
and dark matter, pdm and po do not conserve separately. They must rather enter the energy 
balances [u| 

p D + 3Hp D (l + w D ) = -Q, (32) 
Pdm + 3Hp DM = Q, (33) 
Pbm + 3Hp B M = 0, (34) 

where we have assumed the baryonic matter does not interact with dark energy. Here Q denotes 
the interaction term and we take it as Q = 3b 2 H (pdm + Pd) with b 2 is a coupling constant. 
This expression for the interaction term was first introduced in the study of the suitable coupling 



between a quintessence scalar field and a pressureless cold dark matter field 22j, |23j] . The choice of 
the interaction between both components was meant to get a scaling solution to the coincidence 
problem such that the universe approaches a stationary stage in which the ratio of dark energy 
and dark matter becomes a constant. In the context of HDE models, this form of interaction was 
derived fr o m the choice of Hubble S ca!e as the IR cu t „ ff Q . 

Combining Eqs. (|13j) and (|2ip with the first Friedmann equation ([5]), we can rewrite this 
equation as 

Per + Pk = PBM + PDM + PD + P<j>i (35) 



where we have defined 



P0 = ~aif a 



Dividing Eq. (f35j) by p cr , this equation can be written as 



where 



^BM + ^DM + ^D + ^6 = 1 + &>k, 



Per V 3 



(36) 



(37) 



Thus, we can rewrite the interaction term Q as 



Q = 3b 2 H(p DM + p D ) = 3b 2 H PD (l + r) 



where r = pdm/pd is the ratio of the energy densities of two dark components, 







DM 



-1 + 



1 



i + n k - n BM + 2« ( i - ^) 



Using the continuity equation ([34"]) . it is easy to show that 



3 / 



(38) 



(39) 



(40) 



(41) 



where Qbmo ~ 0.04 is the present value of the fractional energy density of the baryonic matter and 
z = a" 1 — 1 is the red shift parameter. Inserting Eqs. (j23|) . (I39p and (I40D in Eq. (132H we obtain 
the equation of state parameter 

1 2a 2 V / H^ 



WD 



If we define, following 



3c 



cos y — b 2 £lp 



1 + O fc - ft BM + 2a 1 - — 



au; 



(42) 



121 ] . the effective equation of state as 



eff L 
W D =W D + —, 



(43) 



where V = 36 2 (1 + r)H . Then, the continuity equation (|32[) for the dark energy can be written in 
the standard form 



p D +3Hp D (l + wf) = 



(44) 



Substituting Eq. (02]) in Eq. (@3D, we find 



9 



From Eq. (|45|) we see that with the combination of Brans-Dicke field and HDE, the effective 
equation of state, itffj , can cross the phantom divide. For instance, taking Qd = 0.73 for the 
present time and c = 1, the lower bound of Eq. (|45|) is w e p = — ^ — 0.9. Thus for a > 0.15 we 
have w e p < — 1. Therefore, the Brans- Dicke field plays a crucial role in determining the behaviour 
of the effective equation of state of interacting HDE. It is important to note that in standard HDE 
(a = 0) it is impossible to have vf^ crossing — 1 [12j. Returning to the general case (|42p . we see 
that when the interacting HDE is combined with the Brans-Dicke scalar field the transition from 
normal state where wd > — 1 to the phantom regime where wo < -1 for the equation of state 
of interacting dark energy can be more easily achieved for than when resort to the Einstein field 
equations is made. In the absence of the Brans- Dicke field (a = 0), Eq. (|42|) restores its respective 
expression in non-flat standard cosmology [25] 

1 



wd 



3 3c 

Next, we examine the deceleration parameter, q 
Eq. ([28]) . one can easily show 

1 



cos y - b 2 £l D ~ l (1 + Q k - n BM ) . (46) 
-a/(aH 2 ). Substituting wd from Eq. (|42|) in 



2a + 2 

-3b 2 (l + n k 



2 

(2a + l) 2 + 2a(aw - 1) + Q k - (2a + l)tt D - -VL% 2 cos y 



9. 



BAI 



+ 2a 1 



auj 
~3~ 



(47) 



If we take Qd = 0.73 and ~ 0.01 for the present time and c = 1, a ~ 1/w, lo = 10 4 , cosy ~ 1, 
£1 bm ~ 0.04 and b = 0.1, we obtain q = —0.5 which is again compatible with recent observational 
data [34]. When a = 0, Eq. (|47p reduces to the deceleration parameter of the interacting HDE in 



Einstein gravity 



25] 



-(l + n fc )--^-^cos/y 



3/2 
D 



— (1 + VL k — Qbm) ■ 



(48) 



We can also obtain the evolution behavior of the dark energy. Taking the derivative of Eq. ([16 
and using Eq. (|18p and relation (Id = H£l' D , we find 



D 



D 



■ cosy 



(49) 



where the dot is the derivative with respect to time and the prime denotes the derivative with 
respect to x = In a. Using relation q = —1 — we have 



n' D = 2Q D (q + 



cos y 



(50) 
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where q is given by Eq. (|47p . This equation describes the evolution behavior of the interacting 
HDE in Brans-Dicke cosmology framework. In the limit of standard cosmology (a = 0), Eq. (|50|) 



reduces to its respective expression in HDE model 
Q! D = flu 



25] 



I - O l>] j i + M^cosy) — 36 2 (1 - O.um) + £">/. 



(51) 



For flat universe, = 0, and Eq. (|51|) recovers exactly the result of 



2J. 



IV. SUMMARY AND DISCUSSION 

In summary, we studied the interacting holographic model of dark energy in the framework 
of Brans-Dicke cosmology where the HDE density po = 3c 2 /(8ttGL 2 ) is replaced with pn = 
3c 2 (j) 2 /{AuiL 2 ). Here 4> 2 = w/(27rG e fr), where G e fr is the time variable Newtonian constant. In 
the limit of Einstein gravity, G e s — > G. With this replacement in Brans-Dicke theory, we found 
that the accelerated expansion will be more easily achieved for than when the standard HDE is 
considered. We obtained the equation of state and the deceleration parameter of the HDE in a 
non-flat universe enclosed by the event horizon measured on the sphere of the horizon with radius 
L = ar(t). Interestingly enough, we found that the combination of Brans-Dicke and HDE can 
accommodate = -1 crossing for the equation of state of noninteracting dark energy. For 
instance, taking 0_d = 0.73 for the present time and c = 1, the lower bound of wd becomes 
wd = — 0.9. Thus for a > 0.15 we have wd < — 1. This is in contrast to Einstein gravity 
where the equation of state of noninteracting HDE cannot cross the phantom divide wd = — 1 5J]. 
When the interaction between dark energy and dark matter is taken into account, the transition 
from normal state where Wd > — 1 to the phantom regime where < -1 for the equation of 
state of HDE can be more easily accounted for than when resort to the Einstein field equations is 
made. 

In Brans-Dicke theory of interacting HDE, the properties of HDE is determined by the pa- 
rameters c, b and a together. These parameters would be obtained by confronting with cosmic 
observational data. In this work we just restricted our numerical fitting to limited observational 
data. Giving the wide range of cosmological data available, in the future we expect to further 
constrain our model parameter space and test the viability of our model. The issue is now under 
investigation and will be addressed elsewhere. 
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